Accepted for publication in The Astronomical Journal 

Preprint typeset using MTpX style emulateapj v. 12/16/11 



CONSTRAINING STELLAR PROPERTIES OF INTERVENING DAMPED Lya AND Mg II ABSORBING GALAXIES 

TOWARD GRB 050730 

Y. Minowa 1 , K. Okoshi 2 , N. Kobayashi 3 , and H. Takami 4 

Accepted for publication in The Astronomical Journal 



o 



o 

u 

CM 

o 



> 

o 

00 

in 

t> 
o 

<N 



ABSTRACT 

We performed multi-band deep imaging of the field around GRB 050730 to identify the host galaxies of inter- 
vening absorbers, which consist of a damped Lya absorption (DLA) system at z a bs = 3.564, a sub-DLA system 
at z a bs = 3.022, and strong Mg II absorption systems at z a bs = 1 -773 and 2.253. Our observations were performed 
after the gamma-ray burst afterglow had disappeared. Thus, our imaging survey has a higher sensitivity to the 
host galaxies of the intervening absorbers than the normal imaging surveys in the direction of QSOs, for which 
the QSO glare tends to hide the foreground galaxies. In this deep imaging survey, we could not detect any 
unambiguous candidates for the host galaxies of the intervening absorbers. Using the 3-a upper limit of the 
flux in the optical to mid-infrared observing bands, which corresponds to the UV to optical bands in the rest- 
frame of the intervening absorbers, we constrained the star-formation rates and stellar masses of the hosts. We 
estimated the star-formation rates for the intervening absorbers as < 2.5 M yr -1 for z > 3 DLAs and < 1 .0 M 
y r -i f or Zr ^ 2 Mg II systems. Their stellar masses are estimated to be several times 10 9 M or smaller for all 
intervening galaxies. These properties are comparable to dwarf galaxies, rather than the massive star-forming 
galaxies commonly seen in the z > 2 galaxy surveys based on emission-line selection or color selection. 

Keywords: galaxies: evolution - galaxies: ISM - galaxies: stellar content - gamma-ray burst: individual: (GRB 
050730) 



1. INTRODUCTION 

Exchange of gas and metals between the intergalactic 
medium and galaxies is a fundamental process of galaxy for- 
mation and evolution. In the standard model of galaxy for- 
matio n based on the cold dark matter (CDM) c osmology 
(e.g.. IWhite & Reesl [19781: iFall & Efstathioul[l98Qh . galaxies 
are thought to be formed by star formation caused by the 
cooling and condensation of gas in the core of galactic ha- 
los produced by the hierarchical clustering of dark matter. To 
reveal the process of galaxy formation and evolution, it is im- 
portant to explore the gaseous and stellar properties of young 
galaxies at high redshifts. The study of stellar populations 
in high-redshift (z > 2) galaxies has been largely developed 
by the detection of col or- selected galaxie s, such as Lyman 
break galaxies (LBGs; S tei del et all |2003|) . or emission-line 
select ed galaxies, such as Lya emitters (LAEsi lCowie & Hu 
119981) . The study of these galaxy populations provides us with 
the properties of their stellar population, including star forma- 
tion rate (SFR), stellar mass, and age. The gaseous properties 
of high-redshift galaxies, such as hydrogen and metal abun- 
dances, have also been studied in this decade, but these stud- 
ies are mostly based on observations of intergalactic matter 
identified as hydrogen or metal absorption lines in the spectra 
of bright QSOs. It is difficult to connect the stellar properties 
of the emission-line selected galaxies to the gaseous proper- 
ties of the absorption-line selected galaxies, as they have very 
different selection biases. The emission-selected samples are 
biased by their flux and depend on the depth of the imaging 
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surveys. Therefore, they only provide information in galaxies 
at the brightest end of the luminosity function. In contrast, 
the absorption-selected samples are selected homogeneously 
by their gas content, and thus, they are less biased by the flux 
of the host galaxies. To connect stellar and gaseous properties 
of high-redshift galaxies, it is essential to detect the emission 
from the galaxies associated with the high-redshift absorption 
line systems. 

Among the several absorption lines seen in the spectra of 
back ground QSOs, the damped Lya absorption (DLA) sys- 
tems (Wo lfe et alTl 986) present unique opportunities to select 
H i-rich (H I column density (Nui) > 10 203 cm -2 ) galaxies at 
high redshift. The D LAs dominate the H I content of the uni- 
verse at z > 2 (e.g., IPeroux et al.l 12003b iProchaska & Wolfel 
2009: lNoterdaeme et al.ll2009|) and contain a large amount of 
H I gas mass, ac counting for a l arge fraction of the present- 
day stellar mass (ICole et al .11200 1[) . Thus, the high-z galax- 
ies associated with absorbers such as DLAs are expected 
to be progenitors of present-day galaxies. Many attempts 
to identify DLA host galaxies have been made so far. At 
low redshift (z < 1), a fe w dozen DLA absorbing galaxies 
have been established (e.g.. lRao"e t al. 20031: IChen & Lanzettal 
2003; Rao et al. 201 1), although more than 50% of the known 
low-redshift DLAs remain unidentified. At high redshift (z > 
1), however, the success rate for discovering DLA host galax- 
ies is quite low, especially at z > 2. Only a few DLA host 
galaxies have been identified (e.g.. lAMler et al.ll2002l 120041: 
IPeroux et all 12012). This low discovery rate could be caused 
by the fact that most DLA host galaxies are faint and com- 
pact. The glare of the background QSOs also severely ham- 
pers the detection of the faint stellar emission from the DLAs. 
Some independent theoretical predictions suggest that DLA 
host galaxies have typical luminosities of sub-L* with a typ- 
ical proper size of ~ 3 kpc, and 60%-90 % of them could be 
hidden by the glare of brig ht QSOs (e.g.. iFynbo et al.lfl999h 
Okoshi & Nagashima 2005). The limited sample size of iden- 
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tified galaxies associated with DLAs prevents us from contin- 
uing investigations to connect the stellar and gaseous p roper- 
ties of high-redshift galaxies. Recently, Fyn bo et all |2010) 
have initiated an emission-line survey for metal-rich DLAs 
whose luminosities are ex pected to be brighter than typi- 
cal le ss metal-rich DLAs (Mai ler et all 120041: iLedoux et all 
2006). They successfully identifi ed three DLA host galax- 
ies at z > 2 in nebular emission (Fvnbo et al. 2QI3, 120TI 
Noterdae me et al.ll2012l) . As for typical less metal-rich DLA 
host galaxies, only a few have been successfully detected to 
date. Even in the cases of the identified DLA host galax- 
ies, the contamination from the QSO light obscures the de- 
tailed study of stellar properties such as luminosity and mor- 
phology. Some of the identified galaxies have been identi- 
fied only by nebular emission (such as Lya, [O III], Ha). To 
detect the stellar continuum from absorbing galaxies unam- 
biguously and to increase the number of identified hosts, we 
have conducted a deep imaging survey for host galaxies as- 
sociated with intervening absorbers toward gamma ray burst 
(GRB) afterglows. Because the GRB afterglow dims on a 
short timescale, observations in the direction of GRBs sev- 
eral months after their emergence allow for more sensitive 
searches for absorbing galaxies than do those in the direc- 
tion of QSOs, for which the contamination from background 
lights is significant. Many optical spectroscopy observations 
of GRB afterglows in the hours after their bursts have been 
performed to date. Some of these hav e shown signs of in- 
tervening (z a b s < zgrb) absorbers (e.g., [M etzger et al. 1997b 
iMirabal et aLll2002t IVreeswiik et al.ll2004t iChen et al.ll2QQ5h . 
Several attempts have been made to identify the galaxies giv- 
ing rise to the intervening absorption systems in the direction 
of the GRB afterglow. Some of these have successfully iden- 
tified the galaxies associated with DLAs or metal- absorptio n 
systems at z ~ 1 (IVreeswiik et al . 112003k IPollack et al.ll20Q9h . 

In this context, GRB 050730 is a very unique target that has 
intervening DLA features at z > 2, where the discovery rate 
for identifying the DLA host galaxies toward QSOs is quite 
low. GRB 050730 was first detected by the Sw i ft sa tellite 
dHolland et al.ll2005l) on 2005 July 30. IChen et alJ (120051) per- 
formed high-resolution (~ 10 km s -1 ) echelle spectroscopy 
of GRB 050730 4 hr after the initial burst, when the after- 
glow was bright at R ~ 17.7, and confirmed the redshift of its 
host galaxy at zgrb = 3.969 based on a strong DLA feature 
(logA^i ~ 22.15). Additionally, the author s also identified 
four other intervening absorption systems (Prochask a et al.l 
l2007ah : z abs = 3.564 DLA system (logtf(H I) ~ 20.3), 
z a bs = 3.022 sub-DLA system (log#(H I) ~ 19.9), and z a bs = 
1.773,2.253 strong Mg II absorption systems, whose rest- 
frame equivalent widths (W t X2796 ) were ~ 1 A. Chen et al. 
(2005) estimated the metallicities for the intervening DLA 
and sub-DLA systems at z = 3.564 and 3.022 in terms of 
silicon abundance using one of the unsaturated Si II absorp- 
tion lines associated with the each system. They found that 
both of the intervening DLA and sub-DLA systems show low- 
metallicities with [Si/H] < -1.3 and -1.5 ±0.2, respectively , 
which are typical for z ~ 3 DLAs (e.g., iFynbo et all 12008). 
Strong Mg II absorption systems (W r X2196 > 0.3 A) are of- 
ten found to be associated with DLA or Lyman - limit sys- 
tems (logtf(H I) > 17.7; e.g. JChurchiireTalll^OOQtrRao et al] 
2006) and are known to be a strong probe of the intergalac- 
tic medium along the line of sight toward QSOs or GRBs. 
Therefore, deep imaging of the field around GRB 050730 is 
particularly important for detecting or constraining the stellar 



properties of the intervening DLAs or Mg II absorption sys- 
tems at z > 2 for the first time. 

In this paper, we present the results of optical to near- 
infrared deep imaging observations of the field around GRB 
050730. The layout of the paper is as follows: the optical and 
near-infrared observations and data analysis are presented in 
Section 2. In Section 3, we provide the results of our spectral 
energy distribution fitting analysis and discuss the properties 
of the intervening absorbers toward GRB 050730. Our main 
results are summarized in Section 4. In this paper, we adopt a 
ACDM cosmology with ^ A = 0.7, Q M = 0.3 and H = 70 km 
s" 1 Mpc" 1 . 

2. OBSERVATIONS AND RESULTS 

2.1. Optical and Near- infrared Imaging 

We observed the field around GRB 050730 in 2007 May 
and June using the Mul ti-Object InfraRed Camera and Spec- 
trograph (MOIRCS; iSuzuki et al.l [2008) mounted on the 
Cassegrain focus of the Subaru telescope at Mauna Kea, 
Hawaii. We took images in the / and K s bands, with total ex- 
posure times of 3600 and 11160 s, respectively. Seeing sizes 
were roughly 0."5 1 in t he / band and 0."5 8 in the K s band. 
A standard star (FS136; Leggett et al. 2006) was observed for 
flux calibration at the beginning or end of each observation. 
The data were reduced using a purpose-built pi peline software 
package called MCSRED (Tanaka et al. 2011). First, we per- 
formed flat fielding with self-flat frames made from the actual 
science images, and then we subtracted the sky background 
from each image. Finally, the data were co-registered and 
combined. 

We also obtained B-,R C -, and z'-ban d optical images of the 
GRB 050730 field using SuprimeCam (iMivazaki et al.ll2002l) 
mounted at the prime focus of the Subaru telescope. The total 
exposure times for B, R c , and z! bands were 1950, 4000, and 
2530 s, respectively. The seeing size varied between l."2 and 
l."5 throughout the night of the observations. The data were 
reduced with standard procedures (i.e., dark and bias subtrac- 
tion, flat fielding and sky subtraction, distortion correction, 
point- spread function (P SF) matching, and mosaicking) usin g 
the SDFRED package (lYagi et al.ll200l lOuchi etaIJ l2004). 
The images w ere calibrated with observations of the SA104 
standard stars (Landolt 2009) at an airmass similar to the ob- 
servations of the GRB 050730 field. 

2.2. Optical and Mid-infrared Archived Images 

The field around GRB 050730 was imaged with several 
telesc opes to detect its optical afterglow just after the burst 
(e.g.. llPandey et al. 1 120061) . Moreover, subsequent follow-up 
deep imaging surveys to detect its host galaxy were made 
more than several months after the burst w hen the optical 
afterglow has disappeared (Che n et ail 12009). We used the 
archived images to specify the position of the afterglow of 
GRB 050730 and to provide additional leverage to identify 
the candidate of the galaxy associated with intervening ab- 
sorbers. 

We used the 7?-band snapshot images of the field around 
GRB 05073 observed on 2007 July 30 (just after the burst), 
with Very Large Telescope (VLT)/FORS2, obtained from the 
ESO archive. We also used the VLT/FORS2 R and /-band 
deep imaging data observed in 2008 February-May (more 
than six months after the burst; i.e., no afterglow in the im- 
age). The total exposure times of the deep imaging data are 
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Table 1 

Summary of the Observations and Data. 



Band 


Instrument 


Date 


Exp. Time (s) 


Seeing (") 


Limiting Magnitude a 


Subaru imaging 


j 










B 


SuprimeCam 


2007 Jun. 13 


1950 


1.6 


26.21 


V 


FOCAS 


2010 May 5 


180 


0.8 


25.73 


R c 


SuprimeCam 


2007 Jun. 13 


4150 


1.5 


26.48 


z' 


SuprimeCam 


2007 Jun. 13 


2530 


1.5 


25.10 


J 


MOIRCS 


2007 Jun. 8 


3600 


0.5 


25.22 


K s 


MOIRCS 


2007 May 1 


11160 


0.6 


25.04 


Imaging data from archive 










^special 


FORS2/VLT 


2006 Feb.-May 


9800 


0.7 


27.69 


/ 


FORS2/VLT 


2006 Feb. 25 


600 


1.0 


24.96 


L(3.6fim) 


IRAC/Spitzer 


2008, Mar. 11. 


6970 


1.9 


24.65 


Subaru spectroscopy 










300B+SY47 


FOCAS 


2010 May 5 


5400 


0.7 





a All magnitudes are shown as 3a limiting magnitudes measured in 2 x FWHM apertures. The data were corrected for 
Galactic extinction i Schlegel et al. 1998). 



Table 2 

List of the Intervening Absorbers Toward GRB 050730 



Type 


Redshift 


logAf(H I) 
(cm- 2 ) 


Mg II W r X2196 
(A) 


[Si/H] 






L/U a 


SFR 

(M© yr 1 ) 


log(M st ellar/M ) 


DLA 


3.564 


20.03±0.10 




<-1.3 


> 


-20.79 (B) 


<0.22 


< 1.62 b ,2.45 c 


< 9.52 b , < 9.92 c 


sub-DLA 


3.022 


19.90±0.10 




-1.5±0.2 


> - 


-20.49 (V) 


<0.12 


< 1.24 b , 1.87 c 


<9.68 b ,<9.92 c 


Mg II 


2.253 




0.886±0.031 




> 


-19.95 (R) 


<0.08 


<0.75 b , 1.14 c 


<9.66 b ,<9.68 c 


Mg II 


1.773 




0.943±0.023 




> 


-19.49 (/) 


<0.03 


<0.49 b ,0.74 c 


<9.71 b ,<9.54 c 



a The 3cr upper limit in the ^-band image (Tabled, which corresponds to the rest-frame optical at the redshift of the absorbers (B for z - 3.564, V for 
z = 3.022, R for z - 2.253, and / for z - 1.773 ), was used fo r calculating the upper limit of L/L* . The L* magnitudes for B, V, and R bands are drawn 
from a luminosity function at 2 < z < 3.5 from Marchesini et al. 1 2007), and the L* magnitude for the / band is drawn from a luminosity function at z < 2 
from lllbert etall f200l) . 

b The case that no counterpart associated with the absorbers is detected in the deepest R-band image (R > 27.14). 
c The case that the candidate Gl (R ~ 26.69) is associated with the absorber at each redshift. 



9800 and 600 s in the R and / band, respectively. A basic 
calibration including bias subtraction, flat fielding, and sky 
subtraction was performed on each image using the EsoRex 
FORS pipeline. The resultant images have seeing sizes of 

0. "7 and 1 "0 in the R and / band, respectively. 

We also obtained the L-band (3.6/im) Spitzer/IRAC im- 
ages of the field o f GRB 050730 from the Spitzer archive 
( Las kar et aLll201 ll) . The total exposure time was 6970 s. The 
data were reduced using the standard Spitzer Science Center 
pipeline software, whose outputs are Basic Calibrated Data 
(BCD) FITS files for each frame. The image mosaics were 
created from the BCD data using the MOPEX software. The 
FWHM of the PSF in the final image is estimated to be around 

1. "9. 

2.3. Limiting Magnitude 

To estimate the limiting magnitude, we derived the la back- 
ground fluctuation in each band by directly measuring the 
sky fluxes within 2 x FWHM apertures randomly placed on 
the images. Table Q] summarizes the 3 a limiting magnitude 
within a 2 x FWHM diameter aperture and seeing size in each 
band, as well as the observational information. The limiting 
magnitudes listed in this table are corrected for Galactic ex- 
tinction toward the field around GRB 050730 {Ay = 0.156; 
Schlegel et al. 1998). Hereafter, we used the magnitude cor- 
rected for Galactic extinction. 

The 7?-band image of VLT/FORS2 is the deepest among the 
obtained images. The limiting magnitude of the 7?-band im- 
age is 27.69 (AB), which is more than 1 mag deeper than the 



rest. Therefore, we use the 7?-band image for detecting galax- 
ies around GRB 050730. However, we do not use the 7?-band 
image for deriving the luminosity because the 7?-band corre- 
sponds to the rest- frame UV wavelength at the absorbers' red- 
shift (z > 1.7), and thus mainly traces the star formation activ- 
ities in galaxies rather than the stellar properties. The K s band 
corresponds to the rest- frame optical band at the absorbers' 
redshift (B for z abs = 3.969 and 3.564, V for z abs = 3.022, R for 
Zabs = 2.253, and / for z a b S = 1.773), and thus traces the long- 
lived stellar population in the galaxies. Therefore, we used 
the ^-band image to obtain the galaxy luminosity or its up- 
per limit. Note that although the L (3.6/im)-band image also 
traces the rest- frame optical wavelength at the absorbers' red- 
shift, its limiting magnitude (24.65) is brighter than that of the 
K s band (25.04; see Table©. 

We derived the absolute magnitude and luminosity in the 
rest- frame optical for each absorber's redshift corresponding 
to the 3cr limiting magnitude in the K s band (Table [21). The lu- 
minosity is expressed in terms of L* by assuming the L* mag- 
nitude drawn fr om the rest-frame optica l luminosity function 
at 2 < z < 3.5 dMarchesini etal.ll2007L M*(B) = -22.39 for 
Zabs = 3.564, M*(V) = -22.77 for z , hs = 3.022 M*(R) = -22.67 
for z abs = 2.253) and at z < 2 dllbert et al.1 12005L M*(/) = 
-23. 15 for Zabs = 1 .773). We emphasize that the limiting mag- 
nitude for our ^-band image corresponds to a luminosity less 
than 0.3L* at any absorber's redshift, and we note that we 
can reach the sensitivity necessary to detect sub-L* galaxies, 
which a re expected to be major contr ibutors to the DLA cross- 
section ( Okoshi & Nagashima 2005|). 
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Figure 1. Wtf/z'/^^^nO-band images of the field around GRB 050730. The size of each panel is 25" x 25", which corresponds to ~200 kpc x 200 
kpc at the redshift of the intervening absorbers (z a bs ~ 1.5-3.5). North is up, and east is to the left. A red circle in each panel shows the location of the optical 
afterglow of GRB 50730. The radius of the circle is around 10 kpc at the GRB's redshift (z = 3.969), the area in which most GRB host galaxies were discovered 
(Bloom et al. 2002; Prochaska et al. 2007b). No detectable object is seen within 10 kpc around the GRB afterglow, as reported previously in the literature 
(Chen et al. 2009). 



2.4. Sources around GRB 050730 

Figure Q] shows the acquired images of the field around the 
GRB 050730 afterglow. The size of each panel is 25" x 25", 
which roughly covers the area around ~100 kpc from the line 
of sight toward the GRB afterglow at the absorbers' redshift 
fobs ~ 1.7-3.9). We detected 27 galaxies within this field of 
view in the deepest 7?-band image. Only seven out of the 27 
galaxies are detected in the K s band. 

Long-duration GRBs are known to occur exclusively within 
a few kiloparsecs of the center of s tar-forming galaxies 
dBloom et al.ll2002l: iFruchter et al.H2006l) . We confirmed that 
there was no discernible flux associated with the intrinsic 
GRB-DLA absorbers at z ~ 3.9 68 within 10 kpc around the 
afterglow, as already reported in lChen et al.l ([2009), although 
our limitin g magnitude in R band (~ 27.4) is 1 mag deeper 
than that of IChen et al.l (I2009L R c ~ 26.4). 

By contrast, intervening DLA or Mg II absorbers are se- 
lected independently of any emission or stellar population. 
Their sight lines are selected by a cross section of the inter- 
vening galaxies and should preferentially intersect the outer 
regions of the interstellar medium in the galaxies. To iden- 
tify possible candidates for the galaxies associated with the 
intervening absorbers, we first eliminate galaxies that are lo- 
cated too far away from the line of sight toward the GRB af- 
terglow to produce the observed DLA or Mg II absorptions. 
We adopted a scaling relation between the galaxy's luminos- 
ity (L) and halo radius (R) in the form 



R 

R* 



(1) 



to constrain the extension of the gas cloud associated with the 
D LA or Mg II host galaxies following the method described 
in iFynbo et all (119991) . In this formula, the halo radius R* 
and power-law index t are determined by their fit to the ob- 
served DLA or Mg II host galaxies at low redshift (z < 1) as 
R* = 42.86 kpc and t = 0.26 for the DLA (IChen & Lanzettal 
[111, and R* = 107. 14 kpc and t = 0.35 for the Mg II sys- 
tem (IChen et al.ll2010l) Based on the semi-analytic models of 
galaxy evolution (e.g. jMo et al.l l 999). the halo size of galax- 



ies is expected to be smaller at higher redshift, which has also 
been con firmed in several observational studies o f disk galax- 
ies (e.g.. iBarden et al.N2005b iTrujillo et al.ll2006|) . Therefore, 
we used the radius R derived from the above formula as an 
upper limit of the extension of the gas cloud associated with 
DLA or Mg II host galaxies at z > 2.0. Figure [2] shows the ap- 
parent separation between the position of 27 detected galaxies 
in the R band and the position of the afterglow as a function of 
their apparent magnitude (or 3 a upper limit) in the K s band, 
which corresponds to the rest-frame optical band at the red- 
shift of the absorbers. We also plot the upper limits of the 
halo radii associated with the DLA at z = 3.564, sub-DLA 
at z = 3.022, and Mg II system at z = 2.253 and 1.773 based 
on the above formula. The magnitudes corresponding to L* 
are drawn from the rest-frame optical luminosity function as 
in Section 12.31 We found that three galaxies are located be- 
low the limit of the halo radius in Figure [2 We identified 
these three galaxies as candidates for the galaxies associated 
with the intervening DLA or Mg II absorption systems toward 
GRB 050730. The three candidates are marked Gl, G2, and 
G3 in Figures [T]and [2] in order of the apparent distance from 
the position of the afterglow. Gl and G2 are faint and were 
only detected in the deepest 7?-band image, whereas G3 is the 
brightest among the three candidates and was detected in all 
bands. 

2.5. Photometric Color Selection 

To further constrain the candidates for the host galaxies of 
the intervening abs orbers, we used the s tellar p opulation syn- 
thesis models from iBruzual & C hariot ( 2003) and examined 
a possible evolutionary track in the two color (B - R) versus 
(R-K s ) diagrams at each absorber's redshift (Figure[3]). Stel- 
lar ages in the tracks range from 1 Myr to the cosmic age 
at each redshift. Single- stellar population, exponentially de- 
caying star formation with ^-folding time r of 0.05, 0.1, 0.5, 
and 1 Gyr, and constant star-formation models were used in 
the evolutionary tracks. We also used two different me tallici- 
ties, Z and 1/50 Z . We adopted the lSalpeteii (119551) initial 
mass function (IMF) with mass cutoffs of 0.1 and 100 M for 
all cases. No dust attenuation was assumed, as typical DLAs 
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22 23 24 25 

Apparent K g magnitude (AB) 

Figure 2. Apparent K s magnitude vs. separation from the GRB afterglow for 
galaxies detected within the 25 " x 25 " field-of-view around the afterglow. 
Circles show the galaxies detected in the K s band. Squares show the galaxies 
not detected in the K s band but detected in the R band, and then indicate 
the 3a upper limit of the K s -bmd magnitude. The dashed line marks the 
upper limit of the impact parameter (~ 10 kpc) for the host galaxies of GRB 
050730 at z = 3.969. The triple-dot-dashed line, solid line, dot-dashed line, 
and dotted line mark the upper limits of the impact parameter of the DLA 
or Mg II host galaxies at z = 1.773, 2.253, 3.022, and 3.564, respectively. 
We estimated the upper limits of the extension of cold gas clouds associated 
with DLA or Mg II host galaxies using an empirical relation between the 
galaxy luminosities (L) and radius (R), which is described by a power law 
of the form R = 43(L/L*) - 26 (kpc) for D LAs (fChen& Lanzettdl2003h and 
R = 107(L/L*) - 35 (kpc) for Mg II systems dChen et alJ201Ch . 

show extremely low dust attenuation with E(B -V) < 0.04 
(e.g.. lEllison et al. 1 120051: iFrank & Perouxl 12010b . and the at- 
tenuation should be negligible in this two-color diagram. In 
Figured we also show the colors for the candidates, Gl, G2, 
and G3. For Gl and G2, because the exact B and K s mag- 
nitudes were not available, we include the possible ranges of 
color expected from the 7?-band magnitude and the brightest 
limit of the B- and A^-band magnitudes (see Table [I]). We 
found that the color range for Gl (thin gray shaded area in 
Figure [3]) overlaps with any type of evolutionary track at any 
redshift. The color of G3 (circle in Figure [3]) can also be ex- 
plained by the SED expected from the single stellar popula- 
tion tracks at z = 3.022 and z = 3.564. These indicate that Gl 
could be a candidate for the galaxy associated with one of 
the four intervening absorbers, and G3 could be a candidate 
for the galaxy associated with either z = 3.022 (sub-DLA) or 
z = 3.564 (DLA) intervening absorbers. The color range for 
G2 (thick gray-shaded area in Figure© overlaps with the evo- 
lutionary tracks at z = 3.564 and marginally overlaps with sin- 
gle stellar population tracks at z = 3.022 at an age of around 
40 Myr. This also indicates that G2 could be a candidate for 
the galaxy associated with either of the intervening absorbers 
at z = 3.022 (sub-DLA) or 3.564 (DLA). However, the impact 
parameter of G2 is larger than the size limit for z = 3.564 DLA 
and too close to the size limit for z = 3.022 sub-DLA. Thus, 
G2 is unlikely to be a galaxy associated with the intervening 
DLA or Mg II absorption systems toward GRB 050730. 

2.6. Spectroscopy of the Brightest Candidate 

To identify the redshift of the brightest object (G3), we per- 
formed optical spectroscopy with FOCAS mounted on the 
Cassegrain focus of the Subaru telescope ( Kashi kawa et al.l 
l2002h . We used a 0."8 slit, 300B prism, and SY47 order-sort 
filter, which provided a wavelength coverage of 4800-9000 A 



with a spectral resolution of ~ 9.6 A, to identify the rest- frame 
UV spectral features of the absorbers at z = 2.0-3.5, such as 
Lya emission, continuum depression near the Lyman limit, 
and sever al interstellar absor ption features seen in the spectra 
of LBGs (Sha pley et al.ll2003|) . During target acquisition, we 
also took a V-band snapshot image of the field around GRB 
050730, including the three candidates, using FOCAS (see 
Table [T]). The total exposure times were 5400 s and 180 s for 
spectroscopy and snapshot imaging, respectively. Flux cali- 
bration and atmospheric absorption-band correction were per- 
formed using the spectrum of the spectrophotometric standard 
star GDI 53, obtained on the same night. The spectroscopic 
data were reduced in the standard manner using IRAK We 
performed bias subtraction, overscan subtraction, flat field- 
ing with a dome lamp flat, image distortion correction, wave- 
length calibration with a ThAr lamp, background subtraction, 
and sky background subtraction. Finally, the two-dimensional 
spectra were co-registered and combined, and then the one- 
dimensional spectrum of the object was extracted from the 
combined two-dimensional spectrum. In Figure HI we show 
the resultant FOCAS spectrum of G3. We found no obvious 
rest- frame UV features in the obtained spectrum. 

To further constrain the redshift of G3, we also estimated 
the photometric redshift of the galaxy based on the photom- 
etry in B, V, R, I,z', J, K s , and L(3.6/im) bands (points 
in Figure [4]) using th e publicly available software EAZY 
(lBrammeretal.ll2008h . EAZY fits the observed SED of the 
galaxy with a combination of galaxy templates. We used the 
default EAZY template set eazy_vl.0, which consists of five 
templates that span the colors of galaxies in the semi-analytic 
model and an additional template to compensate for the lack 
of young and dusty galaxies in the semi-analytic model. Fig- 
ured shows the best-fit SED template, which is in good agree- 
ment with both the photometric points and the observed spec- 
trum. The redshift of the best-fit SED is z ph ot = 0.519!°;^. 
We also show the resulting x 2 as a function of the photomet- 
ric redshift, which reaches a minimum at the redshift of the 
best-fit template (z p hot = 0.519). Therefore, we conclude that 
G3 is most likely at z ~ 0.5 and is not related to the known 
intervening absorbers. 

3. CONSTRAINING THE PROPERTIES OF INTERVENING 
ABSORBERS 

The above arguments suggest that Gl is most likely to be a 
host galaxy of one of the intervening absorbers toward GRB 
050730, whereas G2 and G3 are unlikely to be related to the 
absorbers. In this section, we constrained the SFR and stellar 
mass for the four intervening absorbers at z = 1.7-3.5 toward 
GRB 050730 by assuming that (1) Gl is a host galaxy of the 
absorber at each redshift or (2) there is no detectable galaxy 
related to the absorbers. Finally, we compared the constrained 
SFR and stellar mass of the intervening absorbers at z = 1.7- 
3.5 with the other high-z (z > 2) galaxy populations. 

3.1. Star Fpormation Rate 

We estimated SFRs of the galaxies associated with the 
intervening DLA or Mg II absorption systems at z = 
1.773,2.253,3.022, and 3.564 from the rest-frame UV con- 
tinuum (R-ban d in the observed fra me). We used the empiri- 
cal formulae of Madau et al. (1998) assuming Sabeter (1955) 
IMF and no reddening. In Table [2 we show the SFR for the 
absorbing galaxy at each redshift in case (1), the case that Gl 
is a host galaxy of the absorbers. We also estimated the SFR 
using the 3-a upper limit of the rest- frame UV continuum in 
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Figure 3. Evolutionary tracks in the two-color (B-R) vs. (R-K s ) diagrams at the redshift of the intervening absorbers from Bruzual & Chariot ( 2003) stellar 
population synthesis models. Stellar ages in the tracks range from 1 Myr to the cosmic age at each redshift. Single-stellar population (circles), exponentially 
decaying star formation with ^-folding time r of 0.05, 0.1, 0.5, and 1 Gyr (crosses, boxes, triangles, and asterisks, respectively) and constant star formation 
(diamonds) models are used in these tracks. No dust attenuation is assumed. The solid and dashed lines show the evolutionary tracks for metallicities of Zq and 
1/50 Zq, respectively. Thin and thick gray shaded areas show the range of color expected from the 3a upper limit of B- and K s -bmd magnitudes (see Table [T} 
and the R-bmd magnitude of the candidate Gl (R ~ 26.69) and G2 (R ~ 25.65), respectively. The thick open circle shows the color of candidate 3. 




Figure 4. Left: photometric measurements of the candidate G3 (circles with error bars) and the best-fit SED (thin line) to the photometry data. The optical 



spectrum obtained by FOCAS is also plotted (thick line). The candidate G3 is most likely a low-redshift galaxy at z r 
a function of the photometric redshift. 



0.5. Right: x vame of the SED fitting as 



case (2), the case that there is no galaxy associated with the 
absorbers. In all cases, we found that the SFR for the DLA 
galaxies at z > 3 w as ~ 2.5 M yr -1 or smaller. At z > 2, 
iPeroux etaLl (1201 2|) reported the upper limit of the SFR for 
typical less metal-rich DLA host galaxies (Z < 0. 1Z ), which 
are similar in metallicity to our DLA (z = 3.564) and sub-DLA 
(z = 3.022) samples, using the upper limit of the Ha fluxes 
derived from the VLT/SINFONI IFU data. Our derived upper 
limit for less metal-rich DLA host galaxies is comparable to, 



or eve n more stringent than, that constrained by lPeroux et al.1 
(120121 - 3.6 M yr" 1 assuming Salpeter IMF). 

In the case of z ~ 2 Mg II system (W r X2196 ~ 1 A), we 
reached an SFR upper limit of 1 .0 M yr -1 , but no host galax- 
ies were clearly detected. By contrast, host galaxies of strong 
Mg II systems (W r X2196 > 0.3 A) at z ~ 1 were detected with a n 
SFR of more than 10 M yr -1 (lLove grove & Simcodl20TTh . 
This diffe rence in the SFR bet ween z ~ 1 and z ~ 2 was also 
argued in iBouche et alJ (120121) for the case of a very strong 
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Figure 5. Observed flux of Gl in the R band (cross, R ~ 26.69) and the 3a flux upper limit in BVIJK S L bands (circles, see Table [T} are plotted as a function 
of the observed wavelength. We fit the model SEDs (solid lines) to the observed flux for the galaxies associated w ith the intervening DLA o r Mg II absorbers 
at z ~ 1.773,2.253,3.022, and 3.564. The model SEDs were drawn from the stellar population synthesis model of Bruzual & Chariot (2003) by considering a 
variety of star-formation histories (SFHs) including single- stellar populations, exponentially declining SFRs with ^-folding timescales from 10 Myr to 5 Gyr, 
and constant star formation, and two different metallicities with Z = Zq and O.O2Z0 . The models were computed for ages from 1 Myr to the cosmic age at each 
redshift. No dust reddening was assumed in the model SEDs. Fluxes of the model SEDs were scaled to fit the 7?-band flux of Gl (see the main text for details). 
Using the scaled flux and mass-luminosity ratio of the model SEDs, we constrained the stellar mass of the intervening DLA or Mg II host galaxies. The inset in 
each figure shows the distribution of the stellar masses derived from the model SEDs for each intervening absorber. 



Mg II system with W r X2196 ~ 2.0 A. They concluded that z ~ 2 
Mg II absorbers reside in smaller halos than z ~ 1 systems. 
Although our sample size is only two, our result for the strong 
Mg II sy stems with W r ^ 796 ~ 1 .0 A supports the scenario pro- 
posed in Bou che et al.l (l2012|) . 

3.2. Stellar Mass 

Using the GRB absorbers, we were able to strictly constrain 
the continuum flux of the intervening absorbers, as the back- 
ground GRB afterglow used for identifying the absorbers had 
completely disappeared when we performed the deep imag- 
ing. Thus, we did not need to consider the ambiguity of the 
noise estimation arising from the background PSF subtrac- 
tion. Moreover, we could eliminate the possibility that the 
galaxy associated with the absorbers was hidden in the un- 
subtracted core of the background PSF, even if there was no 
detectable galaxy around the position of the afterglow. These 
unambiguous constraints on the continuum level provide a 
unique opportunity to constrain the stellar mass of DLA or 
Mg II host galaxies at z > 2 for the first time. 

The upper limits on the stellar masses associated with 
the intervening DLA or Mg II absorption systems at z = 
1.773, 2.253, 3.022, and 3.564 wer e derived by fitting the stel- 
lar population synthesis models of iBruzual & Charlotl ([2003) 
to the upper limit for each object's continuum flux. We used 
the same set of stel lar population synthesis models as were 
used in Section 12.51 We derived the upper limit of the stel- 
lar mass at each redshift for two separate cases: (1) we as- 



sumed that Gl was a host galaxy of the absorber, and (2) we 
assumed that no galaxy was associated with the absorber. Fig- 
ure [5] shows the fitted model SEDs at each redshift of the in- 
tervening absorption systems as a function of the observed 
wavelength, assuming case (1). In the same figure, we also 
plot the flux of Gl in the R band and the 3-a upper limits on 
the fluxes in the other observing bands. In case (1), the model 
SEDs were scaled to fit the 7?-band flux of Gl. In this case, 
we excluded any model SEDs whose scaled fluxes exceeded 
any of the 3-a upper limits for the observed bands other than 
the 7?-band. In case (2), the model SEDs were scaled to avoid 
exceeding any of the 3-a upper limits in the observed bands. 
Table [2] summarizes the results of the stellar mass estimations. 
We found that the stellar mass of the DLA or Mg II absorbers 
at z > 2 is less than a few times 10 9 M , w hich is comparable 
to the mass of dwarf galaxies like LMC (lHarris & Zaritskyl 
l2009h . 

It should be noted that our derived upper limits on the stel- 
lar mass for the less metal-rich DLA at z = 3.564 ([Si/H] 
< -1.3) and sub-DLA at z = 3 mi ([Si/H] = -1.5 ± 0.2) are 
consistent with the stellar masses predicted from the observed 
mass-metallicity relation for z > 2 star-for ming galaxies (e.g., 
iSavaglio et aTll2005l: iMalolino et al.H2008l) . 

3.3. Comparison with Other z > 2 Galaxy Populations 

In this section, we compare the measured (or constrained) 
SFRs and stellar masses for the DLA and Mg II host galax- 
ies at Zabs > 2 against those obtained for other z > 2 galaxy 
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Figure 6. Specific star formation rate (SSFR) as a function of stellar mass 
for DLA host galaxies at z = 3.022 (squares), 3.564 (circles) and Mg II host 
galaxies at z = 1.773 (triangles), 2.253 (stars). In the sam ple plot, we show 
the values for LBGs (Reddy et al. 2006; Erb et al. 2006, pluses) and LAEs 
IGawiser et a l. 2006; Nilsson et al. 2007; ILai et al.ll2008l crosses) at z = 2- 
3. Because the candidate of the host galaxy for each absorber is not clearly 
identified, we plot the SSFRs assuming that Gl is a host galaxy of each ab- 
sorber (filled symbols) or that no galaxy in our image is associated with any 
absorber (open symbols). The solid lines show the SSFRs for constant SFRs 
from 0.1 to 1000 M© yr" 1 . 

popul ations, such as LB Gs (e.g.jReddv et al.ll2Q06l; lErb et al. 
120061) and LA Es (e .g., iGawiser et alT l20Q6: Nil sson et al. 
2007; Lai et al. 2008). To compare the distribution of their 
SFRs and stellar masses, we investigated their specific star 
formation rates (SSFRs). The SSFR is defined as the ratio be- 
tween the SFR and the stellar mass, and it is an indicator of 
the intensity of star formation in the galaxy. In Figure [6l we 
plot the SSFR versus stellar mass for our DLA absorbers at 
Z > 3 and Mg II absorbers at z ~ 2. Given that we do not 
identify a galaxy corresponding to each absorber, we con- 
strained the SSFR for each absorber in two separate cases: 
(1) we assumed that Gl was a host galaxy of the absorber, 
and (2) we assumed that no host galaxy was detected in our 
image. For all cases, we found that our absorbing galaxies 
were located between the low mass end of LBGs and LAEs. 
This implies that our sample absorbers, which consist of less 
metal-rich ([Si/H] < -1.0) DLAs and Mg II systems, should 
not reside in the massive star-forming galaxies (like typical 
LBGs). Indeed, our derived SSFRs for the DLA/Mg II ab- 
sorbers are consistent with the SSFRs of z = 1.7-2.5, less 
metal-rich DLAs, as derived from their abundance pattern us- 
ing star formation histories (SFHs) similar to local irregular or 
dwarf starburst galaxies with weak star formation efficiency 
(iDessauges-Zavadskv et~alll2007l: ICalura et al.ll2009l) . 

4. CONCLUSIONS 

We analyzed deep imaging data of the field around GRB 
050730 in optical to mid-infrared wavelengths to identify the 
host galaxies giving rise to the intervening DLA and Mg II 
absorbers at z = 1.773,2.253,3.022, and 3.564. We used our 
own near-infrared and optical imaging data obtained with the 
Subaru telescope, together with archived optical and mid- 
infrared imaging data obtained with the VLT and Spitzer. We 
examined the color, impact parameter, and photometric red- 



shift (for the brightest galaxy only) for the galaxies detected 
in our multi-color images to constrain the candidates for the 
galaxies associated with the intervening absorbers. We found 
no unambiguous candidate for the host galaxies in our data. 
Using the 3 a limiting flux or the flux of the marginally de- 
tected galaxy (Gl) in the 7?-band image (corresponding to 
the UV wavelength in the rest-frame of the intervening ab- 
sorbers), we placed the following constraints on the upper 
limits of the SFR in the galaxies giving rise to the intervening 
absorbers: 2.5 M yr -1 for the DLA and sub-DLA systems 
and 1.0 Mq yr -1 for the Mg II systems. We used the 3-a up- 
per limit of the flux in each wavelength to constrain the stellar 
mass of the absorbing galaxies to be several times 10 9 M or 
lower. We confirmed that the stellar mass limits for the DLA 
and sub-DLA at z = 3.564 and 3.022 are consistent with the 
stellar masses expected from their metallicity, assuming the 
mass- metallicity relation observed at z > 3 (Maiolino et al.l 
2008). Both the SFR and the stellar mass for the intervening 
absor bers show properties sim ilar to dwarf galaxies like the 
LMC (lHarris & Za ritsky 2009). The intervening absorbers in 
the direction of GRB 050730, which consist of the less metal- 
rich DLAs and Mg II systems with W r X2196 ~ 1.0 A, show 
low SSFRs (< 1.0 Gyr -1 ), which are comparable to LAEs or 
LBGs at the low-mass end of their mass function. This low 
SSFR for our intervening absorbers favors the scenario that 
typical DLAs and Mg II systems with low metallicities are 
likely to reside in dwarf galaxies, rather than in massive star- 
forming galaxies like typical LBGs. 

We thank Kentaro Aoki and Takuji Tsujimoto for helpful 
discussions in the early stage of this work and for careful 
reading of the manuscript. This work is mainly based on data 
collected with the Subaru Telescope, which is operated by Na- 
tional Astronomical Observatory of Japan (NAOJ). Also, this 
work is based in part on archival data obtained from the ESO 
Science Archive Facility and Spitzer Space Telescope, which 
is operated by the Jet Propulsion Laboratory, California Insti- 
tute of Technology, under a contract with NASA. 
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